m 
o 
o 



in 

(N 



o 



The kinetic spherical model in a magnetic field 



^— > 

c« 

a 

C 

o 

(N 
> 

O 
m 
o 

c3 



Matthias Paessens a b and Malte Henkel" 

a Laboratoire de Physique des Materiaux, 1 Universite Henri Poincare Nancy I, 
B.P. 239, F - 54506 Vandceuvre les Nancy Cedex, France 

6 Institut fur Festkorperforschung (Theorie II), Forschungszentrum Jiilich, D - 52425 Jiilich, Germany 

Abstract 

The long-time kinetics of the spherical model in an external magnetic field and below the 
equilibrium critical temperature is studied. The solution of the associated stochastic Langevin 
equation is reduced exactly to a single non-linear Volterra equation. For a sufficiently small external 
field, the kinetics of the magnetization-reversal transition from the metastable to the ground state 
is compared to the ageing behaviour of coarsening systems quenched into the low-temperature 
phase. For an oscillating magnetic field and below the critical temperature, we find evidence for 
the absence of the frequency-dependent dynamic phase transition, which was observed previously 
to occur in Ising-like systems. 



-a 
c 

o 
o 



X 
S3 



PACS numbers: 05.20.-y, 05.10.Gg, 75.10.Hk, 75.60.Ej, 02.30.Rz 



1 Laboratoire associe au CNRS UMR 7556 



1 Introduction 



Non-equilibrium critical phenomena are a subject of intense research activity. A common way to reach 
such a situation is through a rapid change of one of the macroscopic variables which enter into the 
equation of state. For definiteness, consider a simple ferromagnet. It may be brought out of equilibrium, 
for example, by starting initially from a fully disordered state and then quench the system rapidly to a 
temperature below the system's critical temperature T c > 0. The resulting ageing behaviour has been 
in the focus of intensive study, see PjJ [21 El HI El for reviews. Another way to reach a non-equilibrium 
state is to start from a magnetically ordered state below T c and then turn on a magnetic field H 
oriented antiparallel with respect to the magnetic order parameter. Then the system will find itself in 
a metastable stable and a magnetization reversal transition towards the stable ground state will take 
place, see (HI Ej for reviews. 

After a quench to below T c , the system undergoes phase-ordering, that is domains of a time- 
dependent typical size L{t) ~ t x l z form and grow, where z is the dynamical exponent. As a conse- 
quence, a system of infinite size slowly evolves towards an equilibrium state, without ever reaching it. 
This evolution is more fully revealed through the study of two-time quantities, such as the two-time 
autocorrelation function C(t, s) and the autoresponse function R(t, s) 



c(t, S ) = (4>(ms)) , R(t, S ) l){n(t)} 



5h(s) 



(1.1) 

h=0 



where <ft is the order parameter, h the conjugate magnetic field, t is called the observation time and 
s the waiting time. Ageing occurs in the regime when s and r = t — s > are simultaneously much 
larger than any microscopic time scale r m i cro . In many systems, one finds in the ageing regime a scaling 
behaviour, see EJ 

C(t, s) = s- b f c {t/s) , R(t, s) = s- l - a f R {t/s) (1.2) 

where a and b are non-equilibrium exponents. For T < T c , b = while a depends on whether there 
are short-ranged or long-ranged correlations in the equilibrium state. For short-ranged correlations, 
a = 1/z, whereas for long-ranged correlations a = (d — 2 + 7])/z jH]. The scaling functions behave for 
large arguments x = t/s^> 1 asymptotically as 

fc(x) ~ x~ Xc ' z , f R (x) ~ x~ x ^ z (1.3) 

where Ac and Xr are the autocorrelation [HI UH] and autoresponse [TT] exponents, respectively. While 
for a fully disordered initial state, it is traditionally accepted that Ac = Xr, for spatially long-ranged 
initial correlations of the form Ci n i(r) ~ |r|~ d-QI (with a < 0) the relation Ac = A# + a has been 
conjectured [TT] . Furthermore, the rigorous arguments of [12] readily yield Ac > (d + a)/ 2. Very 
recently, different exponents Ac 7^ Xr have also been found in the random sine-Gordon model and in 
addition Ac < d/2 violates the rigorous bound mentioned above [T3] . In addition, and going beyond 
these traditional scaling arguments, it has been proposed recently that the dynamical symmetry group 
of ageing systems might include more general transformations than merely the simple dynamical scaling 
as expressed by eq. (jl.2J) . In particular, there is evidence that the dynamical group of ageing systems 
includes so-called local scale transformations related to conformal transformations in time ^Hj- If that 
is the case, the form of the scaling function 

f R (x)=r x l+a - x ^ z (x-l)- 1 - a (1.4) 

is completely fixed (r is a normalization constant) [Tl] I15j. Going beyond phenomenological tests, at 
least for the case z — 2 it can be shown that, given only the covariance of the response functions under 
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scale and also Galilei transformations, then a Ward identity guarantees the covariance under the full 
group of local scale transformations Furthermore, the causality condition R(t, s) = for t < s also 
follows in a model- independent way [TSj. Tests of Galilei invariance require the consideration of space- 
time-dependent response functions, going beyond the autoresponse function R(t, s) of eq. (jl.2|) . Indeed, 
the phase-ordering kinetics of the Glauber-Ising model has recently been shown to be Galilei-invariant 
in the ageing regime [T7j . 

Another central questions in this context is how to characterize whether/when under the conditions 
just described the system is in thermodynamic equilibrium. It is convenient to consider the fluctuation- 
dissipation ratio [THJ HH] 



At equilibrium, the fluctuation-dissipation theorem states that X(t, s) = 1. The breaking of the 
fluctuation-dissipation theorem has been investigated intensively both theoretically (see e.g. [3J EJ 
miEUE]) and experimentally |22l EHl El . 

Here we are interested in the non-equilibrium behaviour associated with the magnetization-reversal 
transition from the metastable to the equilibrium state. This problem has also received intense attention, 
both experimentally (e.g. |2S1 EHl EE) and theoretically, see [2EJ EHl E01 EU E21 E31 El • What can be 
learned from the study of two-time quantities about this process ? Additional insight may be obtained 
by studying the system's response to a time- dependent, e.g. oscillating, magnetic field and we shall 
study whether there exists a dynamic phase-transition at a finite and non-vanishing value of the period 
P of the field |35| I36| . Surprisingly, there is no such transition in the spherical model, although it is 
known to occur e.g. in the 2D Ising model. 

In order to obtain explicit analytical results, we shall study the the effects of a magnetic field in the 
kinetic mean spherical model, to be defined precisely in section 2. This is one of the very few models 
which can be solved exactly in a great variety of circumstances and has been studied in detail in the 
past, either in the context of continuum field theories [371 EHl EHl HOI EH 132 HH] or else in the form of a 
lattice model [SI E3 EHl 113 SHI EE EJ EH] • It is known that in d < 4 dimensions, the spherical model 
yields results distinct from mean-field theory and therefore permits the study of fluctuation effects. In 
addition, we recall that experimental results of the magnetization reversal [213 EH| are usually described 
in terms of an anisotropic Heisenberg model. Recall that the spherical model shares the following 
equilibrium properties with the 0(3) Heisenberg model and which distinguish it from the often-used 
Ising model: 

• it has a continuous symmetry (0(n) in the n — > oo limit). 

• there is no equilibrium phase-transition in 2D. 

• the equilibrium specific heat exponent a < in 3D. 

These similarities might suggest that qualitatively the kinetics of spherical and the 0(3) Heisenberg 
models should be closer to each other than either is to the kinetics of the Ising model. Still, the spherical 
model should be considered a toy model certainly not meant to be physically realistic. 

This paper is organized as follows. In section 2, the model is defined and the exact solution outlined. 
All physical quantities can be expressed in terms of the time-dependent solution g(t) of a nonlinear 
Volterra integral equation. In section 3, the solution of this equation and its asymptotics are discussed. 
In sections 4 and 5, single- and two-time observables are calculated for the full time-range of the 
magnetization-reversal transition for constant magnetic fields and in section 6 time-dependent fields 




(1.5) 
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are considered. Section 7 presents our conclusions. In the appendices, we comment on the numerical 
techniques and study the exact long-time behaviour of the Volterra equation. 



2 Model and formalism 



We begin by recalling the definition of the kinetic mean spherical model, using the formalism as exposed 
in jlU H3J H7J [TT] . We consider a system of time-dependent classical spin variables S x (t) located on the 
sites cc of a <i-dimensional hypercubic lattice. They may take arbitrary real values subject only to the 
mean spherical constraint 

J] <^(t) 2 ) = Af ( 2J ) 

X 

where M is the number of sites of the lattice. The role of imposing the spherical constraint either 
microscopically or rather in the mean (which is the only case where the dynamics can be solved) has 
been carefully studied recently [13] • Provided the infinite- volume limit is taken before the long-time 
limit, either way of treating the spherical constraint leads to the same results. 

The spherical model Hamiltonian reads 

n = ~ J £ S x (t)S y (t) -J2 H *(t)S x (t) (2.2) 

<x,y> x 

where H x (t) is the space- and time-dependent external magnetic field. The first sum extends over 
nearest-neighbour pairs only and the second sum over the entire lattice. We choose units such that 
J = 1. The system is supposed to be translation-invariant in all directions. The kinetics is assumed to 
be described in terms of a Langevin equation 

= E S v® - ( 2d + i(*))Sm(t) + H x (t) + Vx (t) (2.3) 
y(s=) 

where the sum over y extends over the nearest neighbours of x. The gaussian noise i] x (t) describes that 
the model is in contact with a heat bath. It is characterized by a vanishing ensemble-average and the 
second moment 

(Vx(t)Vy(t')) = 2T5 x J(t-t') (2.4) 
Finally, the function %(t) is fixed by the mean spherical constraint (J2.1|) and has to be determined. 
By a Fourier transformation 

f(q) = T, Ue- iq r , f r = (2n)- d [ dq f(q)e^ r (2.5) 



where the integral is taken over the first Brillouin zone B, the Fourier-transformed spin variable S(q, t) 
becomes 



S(q,t) 



e -u(q)t 



S(q,0)+ I dt'e^'y^if) \H(q,t') +7}(q,t') 
Jo L 



(2.6) 



with the dispersion relation 

d 

( 9 ) = 2£;(l-ooB(ft)) (2.7) 



1=1 
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and we have also denned 



g(t) = exp 2 / dt'$(t 



(2.8) 



Clearly, the time-dependence of S(q, t) and any correlators will be given in terms of the function 

g = g(t). 

We now derive the expressions for the correlators and response functions for an arbitrary external 
field H x (t) and general initial conditions. Consider the two-time spin-spin correlation function 

C x , y (t, s) = C x _ y (t, s) = (S x (t)S y (s)) = (2n)- 2d [ dq dq' ^^-v) (s(q, t)S{q', s)) (2.9) 

Jb 2 

A straightforward calculation gives 
C(q,q';t,s) = (S(q,t)Stf,a)) 

p -w{q)t-u(q')s r / ft \ 

= r---- (27r) d 5 d (q + q')(c(q,t) + 2T dt'e 2 »Wg(t')) 

Vg{t)g{s) L V Jo J 

+(S(q',0)) ftfert&y/rffiH{q,ir) + (S(q,0)) fds' e"«>' ^gjj) H(q', s') 
Jo Jo 

+ f dt' '^/gJf)H(q,t') [* ds' e" (q ' )s ' ^g~{/) H(q', s') (2.10) 
Jo Jo 

where C(q,t) is the single-time correlator. Here the average was carried out over the noise and the 
initial conditions S^O) such that 



(S(q,0)) = J2(S x (0))e' iq - x = (2n) d 6 d (q) S 

X 

In direct space, the two-time autocorrelator becomes 

C(q,q';t,s) 



(2.11) 



C x , x (t,s) = (2n)- 2d / dqdq'e^ +q '> x 

Jb 2 

1 



A[ t -±1)+2T j^duf( t -±^-u)g{u) 



\fg(f)g(s 

+S [ dt'B x (t')^/gJF)+S fds'B x {s')^gJI' 
Jo Jo 

+ f dt'B x {t')^g~{V) [ ds'B x (s')^/gJ^ 
Jo Jo 



(2.12) 



where we have defined 



fit) 


= W d 


fdqe- 2 ^ t =(e- A %(4t)) d 

IB 




A(t) 


= (2ir)~ d 


[ dq e -M9)i C(q, 0) 
IB 


(2.13) 


B x {t) 


= (2vr)- d 


[ dq e^ q)t+iq - x H(q, t) 
IB 





and J is a modified Bessel function [HI]. We see explicitly how the initial magnetization 5*0 and the 
initial correlator C x (0) affect the dynamics of the system. 



It remains to determine the function g(t). Because of the spherical constraint (|2.2j) and spatial 
translation invariance, the equal-time autocorrelator must satisfy 

C (t, t) = [dq C(q, t) = (S x (t) 2 ) = 1 (2.14) 

J B 

This in turn fixes lit) or via (J2.8)) the function git) as the solution of a nonlinear Volterra integral 
equation 

g(t) = A(t) + ^j\t' f{t - t')git') + 2S f\t' B x it')yfg~if) + U\t' B x it')^gJF)^ (2.15) 

For S = and T = 0, eqs. (|2. 1212.1 5 j) had been derived before for the spherical spin-glass [451 • Besides 
on time, git) also depends on the temperature T and the initial conditions parametrized by So and 

O(i). 

The expressions for A(t) and B x it) simplify in certain cases. For uncorrelated initial conditions 

C x j0) = {l-S 2 )5 x , y + S 2 (2.16) 
Then C{q, 0) = 1 - S 2 + (27r)^%)S 2 and 

A(t) = (1 - S 2 ) fit) + S 2 (2.17) 
For a spatially uniform magnetic field H x it) = Hit) we have 

B x it) = Hit) (2.18) 

These two conditions and consequently eqs. (|2.17I2.18|) will be used throughout this paper. 

When So ^ 0, it will be useful to consider besides C(t, s) also the connected two-time autocorrelator 
(see PI] for an analogous situation in the ID Glauber-Ising model) 

T(M) = (S x it)S x is)) - (S x it))(S x is)) (2.19) 

Finally, the response function is obtained in the usual way [SHI EH EJ E51 S3 IHj by considering the 
linear response to the magnetic field. It is easy to see that in Fourier space 



Sh(q,s) 



h r =0 



e -"(*)(*-)J^fl (2.20) 



From these expressions, the autocorrelation function C(t,s) = Co{t,s) and the autoresponse function 
i?(t, s) = Roit, s) can be obtained by integrating over the momentum q. 

Summarizing, the physically interesting correlation and response functions are given by equations 
(|2. 1212. 1912. 20j) together with the constraint eq. ()2.15|) . This constitutes the main result of the general 
formalism. 

In the next section, we turn towards the solution of these equations. Compared to the case without 
an external magnetic field, this task is difficult since the underlying Volterra equation ()2.15j) is nonlinear. 
The mathematical theory of nonlinear Volterra equations is still being developed [52J. In a few cases, 
explicit analytic solutions can be found. Otherwise, we shall turn to numerical methods. 
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3 Solution of the constraint 



It is the peculiar feature of the kinetic spherical model that a complicated many-body problem can 
be exactly reduced to the solution of a single equation. We first derive the exact late-time asymptotic 
behaviour of the solution g(t) of eq. ()2.15|) for a constant field H(t) = H, that is for times t ^> 1. 
Afterwards, we comment on the use of asymptotic expansions for the calculation of physical observables. 

It is convenient to consider first the initial condition So = which is easier to handle. As we shall 
see, the system actually looses its memory of the initial state quite rapidly. 

A first condition on the late-time asymptotics comes from the known fact \C XjX \ < 1. Together with 
eq. (|2.12j) . it is easy to see that the power-law dependence of g{t) on t as found j!7] for the special case 
H = and T < T c is incompatible with that condition in the case at hand. 

We therefore try, for late times t ^> 1, an asymptotic exponential ansatz 

g(t) = ae t ' T , (3.1) 

where a and r are constants to be determined. Indeed eq. (|2.12|) now shows that \C XjX \ is bounded 
if r > 0. To see this, observe that because of the ansatz (|3.1|) the main contribution to the terms in 
eq. (j2.12|) which depend on the magnetic field comes from the upper limit of integration. Consequently, 
the quadratic term in B x {t) dominates over the terms linear in B x (t) and also over those terms which 
do not contain B x [t) at all. For large times t, s we have asymptotically 

lim C x . x (t,s)=AH 2 r 2 , (3.2) 

t,S—rOD 

where the limit is taken for a constant time difference a = t — s > 0. Inserting eq. (|3.1|) into eq. (j2.15|) 
yields for 5*0 = 0, along the same lines 

g{t) = A(t) + 2T fdt' f(t - t')g(t') + AH 2 r 2 g(t). (3.3) 
Jo 

and with f(t) = (e™ 4 * Io(4t)) d from eq. (J2.13|) . Using the Laplace transformation 

POO 

7(p)= / dtf(t)e~ pt (3.4) 

we find from ()3.3j) 

m = JM (3.5) 

This must be consistent with the Laplace-transformed ansatz of eq. ()3.1|) 

W) = ~~~T- (3-6) 

P — T 1 

These two expressions can only be compatible if the denominator in eq. ()3.5|) vanishes at p = r _1 , i.e. 

1 - 2Tj(r- 1 ) = AH 2 t 2 (3.7) 

and this must be a simple intersection (from eq. (j2.17|) we know that A(t) > 0, therefore A(p) > can 
be related to a) . Eq. (J3.7|) is an implicit equation for r and we now show that there is always an unique 
solution, provided H ^ 0. 
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First, we consider the case r —>■ 0. From the definition of f(p) and f(t) > 0, we have f(p) > 0. 
Similarly, f(t) < 1 for t > implies f(p) < p^ 1 . Therefore 

lim (1 - 2Tj(r- 1 )) = 1 (3.8) 

Second, we consider the case r — > oo. From the results of [47J on f(p) one has 

lim (l-2T7(r- 1 )) = 1-^ (3.9) 

It is well-known that the Laplace transformation f(p) of a positive function f(t) decreases monotonously 
with p jJT]. Therefore, the left-hand-side of eq. (|3.7|) decreases monotonously from 1 to 1 — T/T c as 
r increases from to oo, while the right-hand-side ~ r 2 increases monotonously for H ^ 0. This 
establishes the existence of a simple intersection and therefore of a unique r which describes the late- 
time asymptotics of g(t) for H ^ 0, see eq. (|3.1|) . For \H\ — > and T < T c we find r — > oo, while for 
if = a solution for r only exists if T > T c . This reproduces the well-known result that for H = 0, 
<?(£) only has an exponential behaviour for T > T c j57[ [n]. The fact that in the case H ^ we find 
an exponential behaviour for all temperatures T shows that the system relaxes to an equilibrium state 
after the finite time r 03 an d neither critical behaviour nor ageing is expected for late times. It is now 
clear that adding the extra terms coming in for Sq ^ will merely generate sub-leading corrections and 
the asymptotic solution eq. (|3.1j) will not be affected. 

In conclusion, we have established: the leading long-time asymptotic behaviour of g(t) is given by 
eq. where r is the unique solution of eq. \3.7\) and with a = —A(l/r)/(2Tf (1/r)), for any value 

H ^ of the constant magnetic field and any given mean initial magnetization So- 

For finite times, there is no analytical solution of eq. (j2.15|) available. Instead, as described in ap- 
pendix A, we determine g(t) numerically. Although the two-time observables are the relevant quantities 
for the study of ageing phenomena (see section 4), it is still useful to consider single-time observables 
like the average magnetization S(t) given by 



s(t) = (s x (t)) = 1 



git) 



S + fdt'H{t')^ 
Jo 



(3.10) 



In practise, care is required in using asymptotic solutions of g(t) for the prediction of the time- 
dependence of observables such as S(t). We illustrate this in figured where g(t) and the distance of the 
magnetization to its equilibrium value Soo~S(t) are shown as a function of time. We see in figure^ that 
after an initial fall-off, g(t) quickly reaches the asymptotical regime of exponential growth. In figured, 
however, we compare the mean magnetization S(t) as found from the exact numerical solution g(t) (see 
appendix A) with the one obtained from an asymptotic fit of the form 

■Cm ax 

g{t)^e*^a t {t-t )- t (3.11) 
e=o 

where we use £ max = 5. Although that asymptotic fit for g(t) cannot be distinguished from the exact 
numerical result in figure [T^, the deviation in S(t) is considerable. 

In the rest of this paper, we shall use the direct numerical solution of eq. f)2. 15j) . 
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Figure 1: The function g(t) (a) and the distance of the magnetization to its equilibrium value (b) for 
d = 3.5, T = 2 (T c ps 5.27), H = 0.1 and So = 0. The full curve shows the results for the direct 
numerical calculation, the dashed line shows the results for the 5th order fit which coincides with the 
full curve in (a). 

4 Single-time observables 

Our first applications consider single-time observables, which are the ones most commonly studied. 

An instructive example on the importance of fluctuation effects is constructed as follows. For a 
given external magnetic field H, one may easily calculate the equilibrium magnetization M eq . Now 
prepare the system such that the spins have a mean magnetization M eq but such that spins on different 
sites are uncorrelated. The time evolution of S(t) is shown in figure El While a mean-field description 
would have predicted a constant S(t), we see that the magnetization is not constant but increases 
towards a peak before it falls back to the equilibrium value M eq . Intuitively, we would expect that the 
individual spins tend to align with the local magnetic field provided by their neighbours. Since initially 
Sq = M eq > 0, one orientation is preferred with respect to the other one and domains oriented in parallel 
to M eq will grow preferentially. When the domains have grown large enough the influence of this effect 
decreases and the system approaches quickly the equilibrium and the magnetization decreases again. 
This picture, although close in spirit to the Ising model with its discrete spin variables, also works in 
the spherical model, in spite of the fact that the interaction can be reduced to a free-field theory. The 
remnant interaction between different spins provided by the spherical constraint is sufficient to achieve 
non-trivial correlations between different spins. 

We have seen in the previous section that for H ^ and very late times the system relaxes back to 
its unique equilibrium state. For a vanishing magnetic field, the equilibrium free energy would have a 
double-well structure with two equivalent minima, corresponding to the two possible orientations of the 
mean magnetization. Turning on a magnetic field, this potential is tilted and the depth of the two local 
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Figure 2: The magnetization of the system evolving from So = calculated ford = 3.5,T = 2(< T c ), 
H = 0.2 and S « 0.810. 

minima is no longer the same. The lower minimum becomes the unique equilibrium state, the other 
one corresponds to a metastable state. It is clear that if the system is initially prepared in the well 
corresponding to the equilibrium state, it will relax rapidly towards that state. Here we are interested 
how the transition from the metastable state towards the equilibrium state occurs. 

Therefore, we prepare the system with an initial magnetization antiparallel to the given external 
field. In figure we show the time evolution of the mean magnetization. After a short time the system 
reaches the metastable state, independently of the absolute value of the initial magnetization So, and 
where S(t) stays practically constant. The system remains in the metastable state for several decades 
until the magnetization is reversed quite rapidly (although one should not be misled by the logarithmic 
time-scale in this figure which makes the changeover to appear be very fast). In order to understand 
better what is going on we define a characteristic length \(t) of the fluctuations 



where r runs over all sites of the lattice AcZ d and C r (t, s) is the spin-spin correlation function ([2.9)1 . 
The time evolution of A(i) is shown in figure Eb- Starting from a very small initial value, X(t) increases 
towards a maximum value which is reached at the time when S(t) starts to deviate perceptively from 
its value in the metastable state. While S(t) changes its sign, X(t) remains approximately constant at 
its maximal values before it relaxes towards the equilibrium correlation length, with a typical value of 
a few lattice spacings. The coincidence of the times of the reversal of S(t) and the peak in X(t) shows 
that whole domains rather than single uncorrelated spins are flipping. 




(4.1) 
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Figure 3: (a) Average magnetization S{t) for d = 3, T = 2 (< T c ~ 3.96), H = 10~ 3 and S = -0.5 (full 
curve), So = —0.75 (dashed curve) and So = —1 (dash-dotted curve), (b) Squared correlation length 
\(t) 2 of fluctuations for So = —0.5, where the other parameters are as in (a). 

5 Two-time observables 

Having seen that the magnetization reversal passes via an intermediate state with highly correlated 
fluctuations, we discuss in this section how this manifests itself in the behaviour of the two-time quan- 
tities. An important quantity is the time $ after which the magnetization reverses itself. Evidently, 
$ = $(H, T, d), but we have not investigated in detail how $ depends on these parameters in detail. For 
illustration purposes, we shall use in this section the same choice of parameter values as in figure El then 
$ ~ 3000. For finite values of t, g(t) can be readily found from eq. f)2. 15j) using the numerical methods 
described in appendix A. We shall focus on the metastable state by restricting to waiting times s in the 
intermediate time regime s < A magnetization reversal is seen if the initial magnetization is chosen 
antiparallel to the external field. 

Our choice of initially uncorrelated spin with a mean magnetization So can be considered as a special 
case of initially correlated spins. The case of spatial long-range correlations of the form Ci n i(r) ~ |r| -<i-ai 
in the initial state was studied in detail before |3*9*1 ITT] . Formally, this reduces to an initial state with a 
constant mean magnetization in the limit a — > —d. Using the exact results of jTTj for T < T c , we have 

C(t,s) = l-l=M e 2 q 

J c 

R(t,s) = [An(t - s)]- d/2 , (5.1) 

where M eq is the equilibrium magnetization. 

In figure 0^ the correlation function C(t, s) is plotted versus the time difference t — s for several 
values of the waiting times s which are chosen to be in the metastable state, that is s < $ (compare 
figure Hit)- After a short time the curves reach a plateau, with a value very close to the equilibrium 
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Figure 4: The two-time autocorrelation function C(t, s) (a), the correlation of fluctuations T(t,s) 
(b) and the response function R(t, s) (c) plotted vs. the time difference t — s for waiting times s = 
{5; 300; 1000; 2000; 2500}; the data were calculated for d = 3, T = 2 and H = 10~ 3 . In (c) the straight 
line shows the formula [An(t — s)]~ d ^ 2 . 



value C(t, s) = M^ q (a small contribution of the magnetic field can be neglected here). C(t, s) maintains 
itself at this value for approximatively three decades, independently of the waiting time s. When the 
observation time t becomes larger than the magnetization reversal time the correlation function 
C(t,s) changes its sign because the spins at time s before the reversal are anticorrelated to the spins 
at time t after the reversal. However, we point out that the changeover takes more time when the 
waiting time s is increased. The curves rapidly approach the expected value — M^ q because spins of the 
metastable state are compared to the stable state. So we conclude that the correlation function C{t, s) 
is mainly determined by the value of magnetization. 

While C(t,s) measures the time-dependence of the autocorrelation of a given spin, T(t,s), see 
eq. (|2.19|) . measures the fluctuations. This is shown in figure |3Jd. It can be seen that for waiting times 
s < 1000, T(t, s) decreases fairly rapidly as a function of the time difference t — s. In addition, a small 
peak is observed in the region t w But for waiting times closer to the magnetization reversal time $ 
(here for s = 2000 and s = 2500), the fluctuations have become quite substantial and show a larger peak 
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around (t — s) + s ~ d. This may be viewed as another hint for the existence of correlated domains: as 
the spins inside of a domain are highly correlated a fluctuation of a spin within such a domain will cause 
other spins in the domain to follow this fluctuation. In turn, a side-effect of the enhanced correlations 
is a longer lifetime of a spin fluctuation. After the magnetization reversal, T(t, s) rapidly falls to zero. 

Finally, in figure Eh the response function R(t, s) is shown. First, we observe that for a time region 
of at least two decades we recover eq. ()5.1|) . which was derived in [TT] for the case without an external 
field. In this region translation invariance holds and hence no ageing occurs. The system behaves as if 
it were in equilibrium although it is only in a metastable state. Second, for observation times t getting 
closer to the reversal time the response function begins to deviate from this simple behaviour. We 
point out that the curves for all waiting times s still collapse onto each other and that this deviation 
occurs although C(t,s) still has not appreciably changed away from M^. Third, for times t > the 
dependence on the waiting times becomes obvious before the response curves decrease very fast. This 
can be explained by considering that the memory of perturbations is lost during the reversal from the 
metastable to the stable state. 

In order to decide whether the system is in equilibrium or not we shall investigate now the zero- 
field-cooled (ZFC) magnetization which is defined by 

M ZFC = HT I duR{t,u). (5.2) 

J s 

This quantity may be related to the fluctuation-dissipation ratio using eq. (jl.5J) . Because of the 
non-vanishing initial magnetization Sq and the presence of an external magnetic field, the quantities 
C(t, s) and r(i, s) are different and a fluctuation-dissipation ratio is better defined using T(t, s), namely 
X(t,s) = TR(t,s) (dT(t,s)/ds)) \ This had been checked explicitly in the ID Glauber-Ising model 
11 and in certin simple model of glassy behaviour [53] • In spin glasses, it had been shown [T%1 IT^] 
from mean-field theory that X = X(C(t, s)) although that is not necessarily so beyond mean field or in 
simple ferromagnets [HI |!)3J . Nevertheless, this assumption is of good heuristic value. In the spirit 
of the enterprise, let us consider the case where here X = X(T(t, s)). This amounts to saying that T 
serves as a clock for the evolution of the system. Then 

f r(t,t) 

M ZFC /H= / drX(r). (5.3) 

JT(t,s) 

Consequently, when plotting MzFcit, s)/H versus T(t, s) for fixed s (see figure EJ) the slope of the curve 
corresponds to the value of X - provided of course that the assumptions leading to ()5.3|) are valid. 2 
Rather, we find in figure that with increasing waiting time s the curves move from the lower right to 
the upper left. On the other hand, for a given value of s, the system starts in the lower right corner and 
moves rapidly along a curve MzFc(r) — To— T until the metastable value Mzfc/H = 1— M mcta — 1— M eq 
is reached. The slope of unity of this curve is the same as would be found for an equilibrium system. 
Surprisingly, while it undergoes the magnetization reversal, the system then passes through a loop, 
which corresponds to the peak in T(t, s), before it reaches a horizontal line, of height 1 — M eq . The 
movement along the horizontal line is a behaviour typical of the low-temperature phase, indeed through 
the magnetization reversal the system behaves as if the quasiequilibrium branch close to the metastable 
state had to be joined with the low-temperature behaviour after the magnetization reversal. All in all, 
this behaviour is quite analogous to the one observed for Mzfc as a function of C(t,s) for systems 
brought into the two-phase region by a temperature quench |57j . 

2 For metastable systems with detailed balance and for time-scales shorter than the nucleation time, a fluctuation- 
dissipation relation is discussed in |56j . 
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Figure 5: Mzfc/H vs. T, where the parameters are as in figure |U For reference, the grey line gives the 
curve Mzfc/H = IT. In the inset the region of the loops is shown in more detail. 



Of course, all the results in this section depend on having taken s < If we take instead s > 
the system quickly relaxes to its unique equilibrium state. 



6 Dynamics in an oscillating field 

Further aspects of the magnetization reversal transitions become apparent when the response of the 
system to a time-dependent external field H = H(t) is studied. This allows to study hysteresis effects 
- related to the easily measured Barkhausen noise - and has been studied for a long time, see jBl E] 
for reviews. From mean-field descriptions |35J EH1 ISHj; one finds evidence that, depending on the 
amplitude and the period P of H(t), the time-dependent (and periodic) magnetization S(t) = S(t + P) 
changes between two different forms. First, there is a single symmetric solution (corresponding to the 
paramagnetic phase) such that 

S(t + P/2) = -S(t). (6.1) 

Second, there may exist a pair of non-symmetric solutions in the ferromagnetic phase where (|6.1j) does 
not hold. Indeed, the existence of a dynamical phase transition was established beyond mean-field 
theory through simulations in the 2D Ising model with Glauber dynamics [23 EDI EH]- The order 
parameter of this transition is the period-averaged magnetization Q = Q(t) defined as 

Q(t) = jj* dt'S(t'), (6.2) 

where P is the period of H(t). In the Ising model for sufficiently strong fields and/or low frequencies 
Q = and S(t) oscillates around zero, but Q remains finite for smaller fields and higher frequencies and 
S(t) then oscillates around one of the two values of the equilibrium magnetization. Detailed finite-size 
scaling analysis has shown that the exponents of Q(t) and also of the associated susceptibility agree 
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Figure 6: (a) The period-averaged magnetization Q(t) for d = 3, T = 2, Sq = 0, a sinusoidial external 
field with period P = 1.5 and amplitude Hq = 0.2. (b) The Lagrange multiplier g(t). On this scale 
only the behaviour of the bounds can be seen between which oscillations take place; these oscillations 
are shown in the inset. The dashed line shows a power law g(t) = 0.1 t~ 1M m c ■ t~ d ' 2 . 



with those of the equilibrium phase transition of the 2D Ising model [23 HUB EH] • This was further 
backed up by showing that the equation of motion of the order parameter reduces to the 4 -theory with 
noise jHHj (similar studies were also performed on the equation of motion of the anisotropic XY model 

1531 )■ 

Still, this kind of non-equilibrium phase transition need not generically exist. In the g-states Potts 
model with q > 3, for example, a mean-field analysis shows that the time-dependent order parameter 
undergoes a cascade of period-doubling bifurcations, rather than a simple phase transition [20]. It is 
therefore of interest to explore the role of the topology of the phase space further by considering a model 
in a different equilibrium universality class. 



6.1 Behaviour of the magnetization 

We consider the spatially constant but time-dependent external field 

#(t) = tfosin0^Y (6.3) 

The calculation of the observables follows the same lines as in the case of the constant field although a 
larger numerical effort is required. By inserting eq. ()3.1U|) into eq. (j6.2j) the period-averaged magnetiza- 
tion Q(t) is readily obtained. In figure UK a typical example for Q(t) is shown but the behaviour seen in 
this case turns out to be generic. Taking the Ising model as a guide, a heuristic argument |36j suggests 
that a dynamic phase-transition should occur at least for all temperatures and field amplitudes H for 
which a metastable state exists. We therefore used the same values for T as before. However, in the 
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Figure 7: g(t) for d = 3, T = 2 and a sinusoidial external field with Hq = 0.2; the periods are 
P = {0.4; 0.6; 0.8; 0.9; 1; 1.5; 2; 5; 10} (from top to bottom). The dashed line shows a function g mas (t) = 

o.ir 1 - 54 «c.r d / 2 

spherical model we find that for small times Q(t) takes a plateau value before it decays exponentially 
for later times. In principle, and in analogy with the Ising model, one might try to find the dynamic 
phase transition by measuring the time r = r(P, Hq) when the transition between the plateau and the 
decay occurs. Following the practical experience of the Ising model either the scale of H or P can be 
normalized away, see [2H EDI EH] • It should therefore be enough to vary the period P (or the frequency) 
but keep the amplitude H constant and map out t(P,H q ). If there is a dynamic phase transition at 
some critical period P c , the cross-over time should diverge t(P c , Hq) = oo. In practice, however, this 
method is quite slow, because the calculations have to be done for increasingly larger time-scales. 

It is a lot more efficient to study the Lagrange multiplier g(t) which is shown in figure 03- We 
observe that the value of g(t) oscillates between two bounds and the temporal behaviour of the bounds 
correlates with the time-dependence of Q(t). Namely, when Q(t) displays a plateau, the bounds for 
g(t) decay exponentially with time while in the region of the exponential decay of Q(t) the bounds of 
g(t) decay with according to a power law. Therefore, the cross-over time r(P, Hq) can be found by 
determining the intersection of the two regimes for the bounds for g(t). 

In figure [7| we display a typical behaviour of g(t) for several values of the period P. We observe the 
cross-over from a roughly exponential behaviour g exp (t) ~ exp(— t/t) with a relaxation time t towards 
a master curve g mas (t) ~ t~ L54 which is reached for all given values of P for sufficiently long times. In 
principle, one might try to estimate the time of cross-over between these two regimes by looking for the 
intersection of g e xp{t) and g ma s(£) and then further ask when this cross-over time will diverge in order to 
find the critical period P c . Since for finite t this intersection will always occur, a more reliable estimate 
of P c will be given by the condition t -1 (P c , Hq) = 0. 

In figure |S] we show t(P, Hq) for d = 3 and d = 5, that is below and above the upper critical 
dimension of the equilibrium critical behaviour. In all curves, we see that t(P) remains finite for all 
values of P which we considered. Phenomenologically, t ~ 1/P V for P small enough and some exponent 
v > (v ~ 1.4 in 3D; v ~ 1.95 in 5D). The fact that t only diverges as P —>■ is evidence that there 
is no DPT in the spherical model in an oscillating magnetic field, in contrast to established results 
[2*9"1 13*0*1 EH! EHj i n the 2D Ising model and also with results on the n — > oo limit of the 0(n) model 
[BTH IBT] . We also see from figure |H1 that the absence of the DPT is not related to whether or not the 
equilibrium phase transition of the spherical model is in the mean-field regime. 
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Figure 8: t(P, H ) for <i = 3 and d = 5 as read off from g ex p(t), see text. 
6.2 Behaviour of the Lagrange multiplier 

By fitting g m as{t) for d = 3 and d = 5 we find exponents of w = 1.52 ± 0.01 and w = 2.51 ± 0.01 
respectively. From these observations, we conjecture that for sufficiently long times, the Lagrange 
multiplier g(t) satisfies the bounds 

Ci < t w g(t) < C 2 (6.4) 

with an exponent w = d/2 and some constants C\^- Indeed, we have also checked that these bounds 
hold not only for sinusoidal fields H(t), but for triangular and rectangular oscillating fields as well. 
Remarkably, the conjectured exponent w = d/2 of the power-law bounds eq. ()6.4|) coincides with the 
same value found for the kinetic spherical model without a magnetic field [UJ ! We have checked this 
for several values of the dimension d and temperatures T > 0. In appendix B, we derive the bounds 
()6.4j) and especially the exponent w = d/2 in the P — > limit and for T = 0, under mild additional 
conditions. A fully disordered initial state simplifies the calculations but the result remains the same 
for any short-ranged initial correlators. Therefore, the relaxation time t(P, H ) is formally infinite for 
P 1 and T = 0. We have thus shown the absence of a DPT in the physical situation where it would 
have been expected to be seen first. In this respect, the spherical model behaves in quite a different way 
than the Ising model. The rigourous derivation of eq. ()6.4j) is left as an open mathematical problem. 

The absence of a dynamical phase transition is further illustrated in figure El There we compare 
g(t) with a rectangular field H(t) (scaled and shifted for convenience). While for small times, g(t) 
oscillates with the driving period P, we see that with t increasing, an additional peak builds up until 
g(t) oscillates with half the period of the driving field at late times. The fact that g(t) oscillates with 
half the external period P is an indication that the system is described by the symmetric solution, see 
eq. (|6.1|) . The same kind of period-halving is also found for triangular and sinusoidal fields. 

This phenomenon is easily understood: since for small times the magnetization oscillates around a 
non- vanishing value, the global symmetry is broken and the two half-periods of the external field affect 
the system in two qualitatively different ways. However, for later times the magnetization oscillates 
around zero and there is no qualitative difference of the response of the system between the two half- 
periods of the external field any more. This fact is reflected by g(t) actually becoming periodic with 
period P/2, viz. g(t + P/2) = g{t). 

The behaviour of S(t) is illustrated in figures ITUIITT1 and^J For relatively small times (upper panel), 
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Figure 9: g(t) (full line) compared to H(t) (dotted line, scaled and shifted) for different times. These 
calculations were done for d = 5, T = 4, Sq = and a rectangular external field with amplitude 
Hq = 0.6 and with period P — 1. 



S(t) oscillates around the positive equilibrium value and is periodic with period P. It is interesting to 
note that the qualitative shape of S(t) for the rectangular oscillating external field in this regime matches 
closely the one observed in the dynamically ordered phase of the 2D Ising model, see |3U| figure 2(b)]. 
For larger times, the dynamic order parameter Q(t) decreases until S(t) oscillates around zero. The slow 
cross-over towards a solution which satisfies eq. (JbMj) is illustrated in the middle panels of figures HUHT^l 
and in the lowest panels, a situation near to (JbMj) is reached, where S(t) becomes antzperiodic with 
period P/2. In the case of a rectangular field shown in figure El the external magnetic- field amplitude is 
still rather small which results in a linear increase and decrease of the magnetization. For stronger fields 
the magnetization reaches saturation during one half-period and the behaviour of S(t) deviates from 
piecewise linearity. The comparison to the sinusoidal and triangular oscillating external field shows that 
in all three cases the magnetization follows the integrated external field for not too large amplitudes. 

The main result of this section is surprising: in spite of the fact that for T < T c there are just two 
equilibrium states for both the Ising and the spherical models in a (sufficiently small) constant mag- 
netic field, the well-established dynamic phase-transition of the Ising model in a temporally oscillating 
magnetic field is apparently absent in the spherical model. 



7 Conclusions 

In this paper we have investigated the non-equilibrium behaviour of the spherical model in an external 
magnetic field. The model's dynamics is described in terms of a Langevin equation and all quantities 
of physical interest can be expressed exactly in terms of the solution of a non-linear Volterra integral 
equation. In few especially simple cases that Volterra equation can be solved exactly, but we have in 
general used numerical methods. 

First, we studied the magnetization reversal transition, in a temporally constant magnetic field, 
which occurs if the system is initially prepared in near to metastable state from which it relaxes towards 
to unique equilibrium state. We find that the system evolves into the metastable state quickly and 
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Figure 10: Time evolution of the mean magnetization S(t) (full line) compared to the one of a rect- 
angular magnetic field H(t) (dotted line, scaled and shifted) for different time regimes; the parameters 
are as in fig. El 



remains there for considerably long times until it finally relaxes into the stable state. For not too 
large magnetic fields, this transition passes through transient states with long-ranged correlations of 
fluctuations, which means that during the magnetization-reversal transition whole domains rather than 
single uncorrelated spins turn over. 

The two-time autocorrelation function is mainly determined by the magnetization so that connected 
correlation functions, which are more sensible to fluctuations, reveal more information. Again we find 
that the transition involves long-ranged correlations. For times smaller than the transition time i9 
we find an effective equilibrium behaviour although the system is merely in the metastable state. In 
many respects, notably the fluctuation-dissipation relations, we find a close analogy with the ageing 
behaviour encountered in the absence of an external field. But approaching the magnetization-reversal 
the autoresponse function and the fluctuation-dissipation ratio show unusual behaviour, indicating that 
the process is rather complex. Therefore, although the non-vanishing magnetic field H sets a finite 
time scale for the relaxation towards the single equilibrium state, we have found a very rich transient 
behaviour which in many respects is quite analogous to the true ageing behaviour found without an 
external field. 

Second, we looked for a dynamic phase transition in a time-dependent external magnetic field H(t). 
Surprisingly, we find evidence that a dynamic non-equilibrium phase transition, which is known to 
occur e.g. in the Ising model, apparently does not exist in the spherical model. For sufficiently low 
temperatures, we rather find that although the dynamic order parameter Q(t) reaches a plateau value 
for small times, there is always a cross-over to a late-time regime where Q(t) decays away to zero. On 
a technical level, this finding can be represented through the conjecture eq. (16. 4 j) which points to an 
unexpected similarity with the phase-ordering kinetics of the zero-field spherical model. 

Given that several equilibrium properties of the isotropic 0(3) Heisenberg model are closer to the 
ones of the spherical model than they are to the Ising model (see introduction), our results raises the 
question whether a dynamic phase transition for the isotropic 0(3) Heisenberg model in an oscillating 



18 



7.5x10' 
7.0x10 
6.5x10 

101 102 103 104 105 106 

_____ 2x10"' 
M 

-2x10-' 

1501 1502 1503 1504 1505 1506 

2x10"' 


-2x10"' 

4001 4002 4003 4004 4005 4006 

t 

Figure 11: Time evolution of the mean magnetization S(t) (full line) compared to the one of a sinusoidal 
magnetic field H (t) with amplitude H = 0.5 (dotted line, scaled and shifted) for different time regimes; 
the other parameters are as in fig. |_3 

field exists. 3 

Lastly, our results beg the questions what are the effects of a magnetic field on the kinetics of a spin- 
glass and what becomes of the magnetization reversal transition and the dynamical phase transition ? 
However, because of the well-known equivalence [IB] between the spherical spin-glass and the spherical 
ferromagnet, studies in different systems with a true glassy behaviour 4 are needed to shed light on this 
issue. 





3 Existing articles on the DPT in Heisenberg models are either mean-field studies or consider the anisotropic case 
62 (which should be more Ising-like). One might anticipate the existence of a critical n c such that in the 0(n)-model in 
an oscillating field, there is a DPT for n < n c analogously to the Ising model and none for n > n c . 

4 The equilibrium behaviour of the Ising spin glass in a magnetic field has been studied in detail, see |63| and references 
therein. For the spherical spin glass in an oscillating magnetic field short-time numerical calculations give evidence in 
favour of a dynamic phase transition at T = 64 . 
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Figure 12: Time evolution of the mean magnetization S(t) (full line) compared to the one of a triangular 
magnetic field H (t) with amplitude H = 0.5 (dotted line, scaled and shifted) for different time regimes; 
the other parameters are as in fig. El 
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Appendix A. Numerical method 

We briefly discuss the numerical solution of the nonlinear Volterra equation (j2.15j) . adapting standard 
methods [65J to the case at hand. 

Eq. (j2.15|) is cast into the following form, using the equations fl2.17|) and ()2.18|) 

g{t) = (l-S 2 )f(t) + S 2 + 2T jTdf f(t-t')g(t') + 2S j\t' H^^/gJ^+UUt' H{t')^gJ¥]\ . (Al) 

As a first step we will discretize the time by dividing the time interval in N — 1 segments of length k 

t i = ki ; i = 0,1,..., N- 1. (A2) 
The continuous functions f(t) are replaced by the iV dimensional vectors 

/ = (/o,/i,-..,/ J v-l) T , fi = f(ti) (A3) 
and the integrals are replaced by a sum by means of the extended trapezoidal rule jHS] 



x N _ 1 



dxf(x) ~ k 



■t'o 



2^0 + fx + h + • • • + /jv-2 + 2^ N ~ 1 



(A4) 
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Therefore, we have the set of equations 



Fi((g 



F ((g , 
go) , (gi, 



9o),f,H, k) 
gd,f,H,k) 








Fn-i ((go, Vg~o) , (gi, y/gd , • • • , (bn-u V9n-i) , f, H k) = o, (A5) 

depending on the known vectors / and H and the step size k. We have F = g — 1 and 



Fi ((go, y/fto) (gi, y/jji) , f, H, k) = gi 



Tkf + -k'Hf-l 



+ vgl 



SokHt + k 2 [\h ^ + H iV9j\ H i + (! - S o) fi + S o + 2Tk \hi9o + 
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i-l 



i=i 



3=1 



(A6) 



3=1 



3=1 



This set of equations can be solved iteratively: F determines go, F± then leads to g± and so on. However, 
since the Fi are functions of gi and y/gl at each step of iteration two a priori distinct solutions for the gi 
are found. They may be obtained by replacing Gi = ^fgi and solving the resulting quadratic equation 
in Gi. So the question arises which of these solutions has to be used. 

We calculated the two solutions for g when only using either the solution according to the positive 
root ('+'-curve) or the negative one ('-'-curve) — the exact solution should evolve somewhere between 
these two limiting curves. We found that decreasing the step size k results in an approach of the 
curve to the '+'-curve where the latter one only slightly changes. Finally choosing a sufficient small k 
the two curves collapse, so that the exact solution is found. For larger values of k the '+'-curve shows 
only small deviations to the limiting curve, so that in all calculations this solution was used. 

For all calculations a step size of k = 1CT 2 was sufficient, except for the data shown in Fig. 121 where 
k = 10 -4 was used because there the time scale is much smaller. 

The evaluation of the one- or two-time observables proceeds by a straightforward implementation 
of their defining integrals by the extended trapezoidal rule. 



Appendix B. 

We derive the bounds eq. ()6.4|) for the Lagrange multiplier g(t) in an oscillating magnetic field H(t), 
for the special case of vanishing temperature T = 0. For convenience, we consider a fully disordered 
initial state and vanishing initial magnetization So = 0, but our results also hold true for arbitrarily 
short-ranged initial conditions and So ^ 0. We define G(t) := \J g(t) . The non-linear Volterra integral 
equation then is 

G(tf = A(t) + Qf ds H{s)G{s)J (Bl) 

We shall assume throughout that G(t) is bounded on the positive real axis, that is \G(t)\ < M < oo 
for t G [0, oo). This assumption is made plausible by our numerical results displayed in figure 01 
Furthermore, we shall assume that the magnetic field itself is bounded, \H(t)\ < Hq. It then follows 
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from (jBl|) that |G(i) 2 | < A(t) + H^M 2 ^. For long times, we may therefore expect a leading power-law 
behaviour which we may write as G(t) ~ t~ w / 2 . We wish to estimate w (the above argument gives 
w > -2). 

Proposition: Let G(t) be a solution of eq. \B1\) and assume that there is a constant M < oo such that 
\G(t)\ < M for all times t G [0, oo). Furthermore, the oscillating field H(t) is assumed to be bounded 
\H(t)\ < H , piecewise continuous and to have the Fourier expansion 

H(t) = H Y,b n sin -p (B2) 

n=l ^ ' 

such that 

00 _ 

B:=Y^-\b n \ (B3) 
n=l 

zs convergent. Finally, let A(t) = fit) = e~ Mt I (4:t) d , where I is a modified Bessel function. Then the 
exponent w of the asymptotic form G(t) ~ t~ w l 2 for t — ► 00 is for P <C 1 

w = ~ (B4) 



Proof: First, for any mag netic field H(t), we trivially have from eq. (jEJ) that G(t) 2 > = /(*) 
(e _4 *Jo(4t)) d and therefore, as t — > 00 



G(t)> 



??r )d/4 



(B5) 



Consequently, w < d/2. 

Second, we wish to find a sharp lower bound on to. This requires some preparations, however. We 
begin by a discussion of the continuity of G(t). Let e > and consider 



G(t + ef - G{tf = (G(t + e) - G(t)) (G(t + e) + G(t)) 

/i — 1 - £ / /*t~\~s fij 

dsH(s)G(s)(j dsH(s)G(s) + J dsH(s)G(s 



(B6) 



Because of the first part and since A(t) decreases monotonically with t, we have G(t + e) + G(t) > 
2y/ A(t + e) and obtain the estimate 



\G{t + e) - G(t)\ < e 



< e 



\Mt A )\ 



2^/A(t + e) 
\A(t)\ 



+ 



t+e 



ds 



\H(s)G(s)\ 
2y/A(t + e) 



ft+E 

2 ds \H(s)G(s)\+0(e] 
Jo 



2^A(i) y/Mfi 



(B7) 



Here the mean value theorem was applied to A(t) where is some intermediate value, t& G [t, t + e]. 
Taking the limit e — ► 0, we conclude that G(t) is Lipschitz-continuous. Then we can apply the mean- 
value theorem to eq. ()B6j) . Because of the continuity of G(t), the limit 



lim 



G{t + e)-G(t) _ A(t) 
i ~ 2G{t) 



+ H(t) I dsH{s)G{s) 
Jo 



(B8) 
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exists for all times t G [0, oo). Taking the derivative of (jBl|) with respect to t, G(t) satisfies the 
differential equation 

G(t) = ^ + H(t) y/G{tf - A(t) (B9) 

Therefore, if H(t) is continuous, G(t) is also continuous. However, if H(t) has jumps, there may be 
jumps in G(t) as well. 

We shall take the average of eq. ()B9j) over the period interval [t, t + P] of the external field. In order 
to prepare this, let <p(t) be a continuously differentiable function. Then 



P 



t+p 



ds sin 



2im 
~P 



s ip(s) 



P (2jn^ rtt)-*t + P) 1 



27m 
P 
2nn 



P 



2im 



t+p 



ds cos 



2nn 
~P 



2im \ . / 2im , 

cos I -^t (p[ti) + cos — =-*2 <p{t 2 ) 



p 



p 



s ) <p(s) 
(BIO) 



where the mean- value theorems were applied and ti, t 2 are intermediate values from the interval [t, t+P] 
For P —>■ 0, we expect ti,t 2 — > t. Since — 1 < cosx < 1, we obtain the following bound for P 1 



P 



t+p 



ds sin 



2im 
~P 



s ip{s) 



<^(r)\+o(P) 



(Bll) 



Before carrying out the average over eq. ()B9|) . we consider the approximation of H(t) as given by 
eq. ()B2j) through a finite Fourier-sum H^t) := Z~2n=i sin (^t). For every finite value of N, H N (t) 
is continuous and if we use Hjf(t) in eq. (|B9j) . so is G(t). Furthermore, if -ff(t) is continuous, then 
i7jv(i) converges uniformly and the N — > oo limit and the integral may be interchanged, viz. 



\im I dtH N (t)ip(t) = \dtH{t)ip{t) 



(B12) 



where ip(t) is some suitable function. If on the other hand H(t) is only piecewise continuous, we only 
have point-wise convergence H^(t) — > H(t). In this case, the well-known Gibb's phenomenon occurs 
which states that close to jump continuities the trigonometric approximations H^(t) overshoot the limit 
H (t) by about 9% of the jump height, see Because H(t) is bounded, we certainly have the bound 
\HN(t)\ < 3H for all iV e N sufficiently large. Then the conditions of Lebesgue's theorem, see [66J, are 
satisfied and one arrives again at (|B12|) . 



Averaging eq. (jB9|) over a period interval, we obtain 



1 

P 



t+p 



dsG{s) 



1 

P 



H lim 

2G(s) P ^v-oo 



t+p 



dsH N (s)^/G(s) 2 -A(s) 



(B13) 



We now consider Pel. The left-hand side and the first term on the right-hand side are estimated by 
the mean-value theorem. For the second term, we use the inequality (|B11|) term by term, taking the 
N — ► oo limit at the end. This gives, up to terms of order o(P) 



G(r) 



1 

< - 
~ 2 



1 

< - 
~ 2 



Mr) 



G(r) 
Mr) 



G{t) 



+ 



+ 



H PB 



71 
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Mr) 



H Py\b n \ 

n=l 



HnPB 



2y/G(T)* - A(T) 



G(t) 



71 



A(t) 



G(t) (B14) 
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where the condition (jB3|) was used. We finally arrive at 



< 



Mr) 



2G(t) 



H PB A{t) 



71 



G{ty 



where <&(a,x) 



1 + 



l\-x 



(B15) 



We now assume that the strict inequality w < d/2 holds. For large values of t, we then have 
A(t)/G(t) 2 ~ fw~d/2 _^ q Q n foe Q^gj- hand, the function x) has a pole at x c = 1 — a 2 = 
1 — (HqPB/tt) 2 . There is a to sufficiently large, such that, say, A(t ) /G(to) 2 < x c /2 and for r > to, 
we have Q(a, A(t)/G(t) 2 ) < x c /2), which is finite. In addition, for uncorrelated initial conditions 
A{t) = f(t), thus A{t) = Ade- Mt Io(At) d ~ l \h(Adt) - I (4dt)] ~ r rf/2_1 . Inserting these asymptotic 
forms into eq. (|B15|) . we find d/2 < w in contradiction with the assumption w < d/2. q.e.d. 

The condition (JB3)) is trivially satisfied for a sinusoidal field. For a triangular field, one has b n ~ n~ 2 
and for a rectangular field, b n ~ In both cases B is a finite constant and we have w = d/2, in 
agreement with our numerical observation. 

The condition P 1 is essential. For P finite, numerical calculations show that G(t) goes towards 
a constant, modulated by a periodic function, thus w = 0. As P becomes smaller, an intermediate 
regime appears, where G(t) ~ t~ d ^, up to a modulation, before the regime mentioned above is reached. 
This will be described in detail elsewhere. 

For any short-ranged initial condition, A(t) ~ t~ d l 2 and our result w = d/2 stays the same. On the 



other hand, for long-range initial correlations of the form C in \{r) ~ |r 
in the same way w = (d + a)/2. 

Finally, for a non- vanishing initial magnetization So ^ 0, we have 



-d—a 



with a < 11, we obtain 



G(t) 



Mt) 

2G{t) 



+ SoH{t) + H(t)^G(t) 2 -A(t) 



(B16) 



Since f* +P ds H(s) = 0, the exponent w = d/2 is unchanged. 
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